Two-dimensional trapping of dipolar molecules in time-varying electric fields 
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Simultaneous two-dimensional trapping of neutral dipolar molecules in low- and high-field seeking 
states is analyzed. A trapping potential of the order of 20 mK can be produced for molecules like 
ND3 with time-dependent electric fields. The analysis is in agreement with an experiment where 
slow molecules with longitudinal velocities of the order of 20 m/s are guided between four 50 cm long 
rods driven by an alternating electric potential at a frequency of a few kHz. 

PACS numbers: 33.80.Ps, 33.55.Be, 39.10.+j 



Cold molecules offer new perspectives, e.g. for high 
precision measurements [lj and collisional physics stud- 
ies 0- Pioneering work on cold molecules has been done 
using cryogenic buffer gas cooling [3( . Another promising 
technique for the production of cold molecules is based on 
the interaction of dipolar molecules with inhomogeneous 
electric fields. For example, low-field seeking molecules 
(LFS) have been slowed down in suitably tailored time- 
varying electric fields Q and have been trapped in in- 
homogeneous electrostatic fields jjj 0- Furthermore, 
efficient filtering Q of slow LFS from an effusive thermal 
source using a bent electrostatic quadrupole guide has 
been demonstrated 

Compared to LFS, the manipulation of high- field seek- 
ing molecules (HFS) is much more difficult. This is 
mainly due to the fact that electrostatic maxima are not 
allowed in free space, and hence, HFS are quickly lost 
on the electrodes. Nevertheless, guiding of HFS in Kep- 
pler orbits |?j and deceleration as well as acceleration 
of HFS ^(j is possible. Despite this progress in manip- 
ulating dipolar molecules, all techniques realized so far 
are suited only for either LFS or HFS, not both simul- 
taneously. However, in future experiments with trapped 
samples of cold molecules, collisions or the interaction 
with light fields are likely to change HFS into LFS and 
vice versa. Therefore, a technique to trap both species si- 
multaneously is vital. In this Letter, we investigate both 
theoretically and experimentally a new technique, which 
can trap both HFS and LFS. In particular, we report on 
the first experimental demonstration of two-dimensional 
trapping of slow ND3 molecules from an effusive source 
in a bent four- wire guide driven by an alternating electric 
field. 

Trapping neutral particles in oscillating electric fields 
works as follows ll|. The force on a molecule in an 
electric field, E, is given by F = -VW(E), with W(E) 
the Stark energy of the molecule. Polar molecules like 
ND3 and H2CO predominantly experience a linear Stark 
shift, W — s E, where s is the slope of the Stark shift. 
Other molecules, however, and atoms usually experience 
a quadratic Stark effect, W = — \aE 2 , with the polar- 
izability a. The inhomogeneous electric fields required 
for trapping can be realized in many different ways. Our 
trapping configuration is sketched in Fig. [I] It consists 
of 4 parallel rods with voltages as indicated. The field is 
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FIG. 1: Schematic of the four- wire setup for typical operating 
voltages of ±5kV, where E — 35kV/cm and f3 — 2.3 x 
10 3 kV/cm 3 . The voltages on a pair of opposite electrodes are 
switched as shown in (a) in order to produce a time-varying 
field that alternates between configurations 1) and 2) shown 
in (b) with a repetition rate in the kHz range. The radius 
Vmax used in simulations is indicated by the dashed circles. 



rapidly switched between two dipole-like configurations 
with angular frequency oj. Close to the center, the field 
can be expanded harmonically, 



E = E + H(t)f3(x 2 -y 2 ), 



(1) 



where Eo is the field in the center and (3 is (half) the 
curvature of the field. The step function H(t) = 1 if 
nT < t < (n + \)T and H(t) = -1 otherwise, with 
T = 2-k/uo the period of the driving field and n an integer. 

Independent of the slope and the sign of the Stark shift, 
the saddle-like dipole potential derived from Eq.(l) con- 
fines the particle in one direction and repels it in the 
perpendicular direction, depending on time. Therefore, 
the time average of the force is small at every position, 
and identical to zero for a linear Stark shift. However, 
the particle performs a micromotion which is locked to 
the external driving field so that the time-averaged force 
does not cancel and the particle experiences a net at- 
tractive force towards the center. This conclusion is in- 
dependent from the exact form of H(t). Therefore, a 
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sinusoidal change between two configurations would also 
work. However, instantaneous switching is more conve- 
nient to realize in the laboratory. 

It is not difficult to solve the equation of motion for 
molecules in this dynamic field. As the x and y degrees 
of freedom separate, the equation of motion is reduced 
to one dimension. It reads x = —H(t)Q 2 x, where we 
have introduced Q 2 = 2s (3 /m for a linear Stark shift 
and Vt 2 = 2aEo/3/m for a quadratic Stark shift of a 
particle with mass m. Following pjj, we search for 
orbits of the form (x(t),x(t)) = U(t)(x(0),x(0)) with 
U(nT + t) = U(t)[U(T)] n . Periodic solutions have the 
form U(T) = U 2 (T/2)U 1 (T/2), with (U^t), U 2 (t)) = 
ff cosh(tft) ^sinh(£f})\ / cos(tfi) ^ sin(tft) \\ 
^fisinh(tfi) cosh(tfi) ^ ^-fisin(tfi) cos(tfi) J J 
the evolution operators for the first and the second half 
cycle. From this, the main stable region can be found 
for |f2T| < 3.75. Narrow higher-order stable regions also 
exist but will further be neglected. Hence there exists 
a sharp drive- frequency threshold for every Stark shift, 
above which trapping is possible, irrespective of the sign 
of the Stark shift and independent of the initial condi- 
tions. 

Note that the finite extension of our trap, \x\ < Xmax, 
limits the trap depth. This can be estimated by averaging 
over the micromotion, which is assumed to be smaller and 
faster than the macromotion. This results in a trap depth 
m7r 2 ^ ax l] 4 /(24w 2 ), typically of the order of 10 mK Q. 
For a more realistic estimate, a detailed analysis is made 
below. 

It is interesting to compare molecule trapping with ion 
trapping in a Paul trap 14]. An ion experiences the 
Coulomb force Fq, which has zero divergence, V -Fq = 0. 
As a consequence, the the magnitude of the force is the 
same for positive and negative ions. For molecules with 
a linear Stark shift and in the harmonic approximation, 
V • F = 0, and there is an accidental correspondence 
with the equations of motion of an ion in a Paul trap. 
It should nevertheless be emphasized that trapping neu- 
tral molecules differs from ion trapping in many ways. 
First, whereas ions are monopolar particles which can be 
trapped in an oscillating quadrupole field, neutrals have 
permanent or induced dipoles which can be trapped in 
an oscillating dipole field. Second, in contrast to ions, 
the force on a molecule in an electric field depends dra- 
matically on the internal state of the molecule. Third, 
as is obvious from Fig. ^ the harmonic approximation 
is only valid in a small region around the center. In 
contrast to ions, which have vanishingly little kinetic en- 
ergy compared to the 10 5 K trap depth, the temperature 
of the molecules is likely to be comparable to the trap 
depth. Hence, molecules explore the full volume of the 
shallow trap, including the outer regions. These outer 
regions are of paramount importance in our experiment. 
It is here that the harmonic approximation breaks down. 
The molecule-field interaction even allows V • F ^ 0. 
This leads to different behavior for HFS and LFS away 
from the center, where the time- average of the force in 



configurations 1) and 2) is not zero. In particular, the 
regions near the electrodes are net repulsive for LFS and 
net attractive for HFS. 

To take all these effects into account, we have per- 
formed a two-dimensional Monte Carlo simulation. In 
the simulation, point-like particles are injected on-axis 
in a random phase of the driving field. The input veloc- 
ities v x and v y are varied. The particles are propagated 
under the influence of the periodically poled Stark force 
for 10 ms, a typical time in the experiment. They are con- 
sidered lost if they leave a radius r max = 1.25 mm. The 
amount of trapped trajectories can be expressed as an 
area A in the two-dimensional v x ^-velocity space quad- 
rant. A is a good measure for the flux, if the guidable 
velocities are equally present in the source. Results for 
four different molecular states of ND3 exhibiting linear 
Stark shifts with a slope, s, of ± 0.6 cm _1 /(100 ^) and 
±1.2^-7(100^) are plotted in Fig. For a molec- 
ular state with a Stark shift of 1.2 cm _1 /(100 — ) and 

/ v cm ' 

at 9 kHz for our experimental parameters, the occupied 
1^2/ -space has almost a quarter-disc shape, with which 
a trap depth of approximately 20 mK can be associated. 
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FIG. 2: Area of guided molecules in transverse- velocity space 
obtained from a two-dimensional Monte- Carlo simulation. 
Curves for LFS (▼, ▼) and HFS (A, A) ND 3 molecules with 
the given Stark shifts are displayed. 



As expected from the analytic analysis, there is a sharp 
turn-on frequency above which guiding is possible both 
for LFS and HFS. For increasing frequencies, the area 
A reaches a maximum before decreasing again for higher 
frequencies. Note that the maximum of A is approx- 
imately four times larger for LFS than for HFS, indi- 
cating that LFS can be guided more efficiently. This is 
due to the anharmonicity of the trapping potential away 
from the center, where stable trajectories exist for LFS 
with comparatively high initial velocities. This effect also 
causes the substructure in the LFS curves in Fig. For 
even higher frequencies, the amplitude of the micromo- 
tion becomes very small and the force on the molecules 
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approaches the time-average of the force in configurations 
1) and 2). This time- averaged force is small near the 
center but substantial and repulsive near the electrodes. 
Stable orbits therefore exist in the high-frequency limit 
for LFS. The stable orbits occupy disjunct areas in veloc- 
ity space, of which the area, A, saturates. The notion of 
a trap depth becomes difficult in this limit. The behavior 
of HFS, however, is completely different: With the dis- 
appearing micromotion, there are no stable trajectories 
and A vanishes in the limit of high drive frequency. 

Obviously, the motion in the full anharmonic poten- 
tial can be very complicated. In fact, we have found nu- 
merical evidence for chaotic motion for some initial con- 
ditions. But despite of the complicated dynamics, LFS 
and HFS are simultaneously trapped and there is a broad 
overlap between the trapping regions for different Stark 
shifts at the same frequencies, as evident from Fig. [21 
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FIG. 3: Schematic of the experimental setup. ND3 molecules 
emerging from the nozzle are injected into the guide. The 
slowest molecules are kept within the guide and after passing 
two differential pumping stages they enter a UHV chamber 
where they are detected with a mass spectrometer. The fast 
molecules are pumped away. 



Our experimental demonstration of AC guiding uses 
an effusive source of thermal polar ND3 molecules, see 
the setup shown in Fig. [3| 8]. The double-bent guide 
consists of four wires which pass through two vacuum 
chambers with two differential pumping sections before 
ending in a UHV detection chamber. The guide has a 
length of 500 mm and is made of 2 mm diameter stainless 
steel rods, with a 1 mm gap between neighboring rods. 
The radii of curvature of the two bends are 25 mm and 
the rods are built around the 0.8 mm inner diameter ce- 
ramic nozzle, constituting the effusive source. Typical 
operation pressures in the nozzle are around 0.05 mbar 
in order to maintain molecular-flow conditions. Most of 
the molecules are not guided and escape into the first 
vacuum chamber, where an operational pressure of a few 
times 10 -7 mbar is maintained. In the detection cham- 
ber, where a pressure below 10 -10 mbar is achieved, the 
guided molecules are detected with an efficiency of about 
10 -4 counts /molecule by a quadrupole mass spectrome- 
ter (QMS)[Hiden Analytical, HAL 301/3F] with a pulse- 



counting channeltron. The pulses from the channeltron 
are processed in a multi-channel scaler. The ionization 
volume of the QMS begins 22 mm behind the guide. To 
protect the QMS from the high electric fields, a metal 
shield with a 5 mm diameter hole is placed 1 cm behind 
the exit of the guide. The time- varying electric fields are 
generated by switching high voltages with fast push-pull 
switches. Switching frequencies are typically in the range 
of a few kHz. In each phase of the alternating field, op- 
posite rods carry voltages of up to ± 7 kV. Due to the 
inversion splitting of the vibrational ground state, ND3 
shows in good approximation a linear Stark shift in the 
relevant electric field range (0-100) kV/cm [TEj . 
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FIG. 4: (a) Molecular flux as a function of the applied fre- 
quency for rod voltages of ±5kV. The graph in the inset 
shows the frequency of maximum flux, cjo, as a function of 
the applied voltage. The dashed line indicates the ujo oc y/V 
dependence, (b) Flux as a function of the applied voltage. 
For each voltage, the frequency is adjusted for maximum flux 
according to (a). The line is a quadratic fit. 



For detecting ND3 molecules the QMS is set to mass 
20. The influence of other gases at this mass is negli- 
gible. The QMS does not provide state-selective detec- 
tion, which is very difficult to achieve if one considers the 
variety of molecular states involved in a thermal (room- 
temperature) ensemble. Our simulation, however, shows 
that the detection of HFS is suppressed, because under 
the influence of the high electric fields most of the HFS 
turn around at the exit of the guide, moving backwards 
on its outside, or they are reflected back into the guide. 

To discriminate the guided molecules from background 



gas and dark counts, the alternating electric field is peri- 
odically applied for 200 ms followed by a 200 ms where no 
field is applied. This on-off sequence is repeated about 
40000 times to obtain a good signal-to-noise ratio. From 
the instant when the alternating field is switched on, 
the molecules can propagate along the guide and the 
signal rises. After about 25 ms, the signal has reached 
50% of the maximum increase, corresponding to veloci- 
ties around 20m/s or about 0.5 K. When the alternating 
field is switched off the flux instantaneously falls off. The 
flux of guided molecules is determined by subtracting the 
background count rate from that during the last 50 ms of 
the 200 ms- long interval where the field is applied. The 
same measurement without injecting gas yields no signal. 

We now investigate the frequency dependence of the 
system. For electrode voltages of ±5kV, the switching 
frequency is varied between (4-15) kHz, and the flux of 
guided molecules is recorded in frequency steps of 500 Hz. 
The experimental data in Fig. 2t show that the signal 
amplitude reaches a pronounced maximum at a frequency 
of ~8 kHz. No discernible signal is observed for fre- 
quencies below 4 kHz. This behavior can be understood 
from the simple analytic theory above, acknowledging the 
thermal ensemble of molecular states. For every molec- 
ular state there exists a cut-off frequency correspond- 
ing to the associated Stark shift below which no guid- 
ing is possible. This frequency increases with the Stark 
shift. For example, the cut-off frequencies for molecules 
with Stark shifts of 0.6 cm -1 and 1.2 cm -1 at an electric 
field strength of lOOkV/cm in our guide are 3.3 kHz and 
4.7 kHz, respectively. For small frequencies stable tra- 
jectories exist only for molecules with small Stark shifts. 
These molecules are unlikely to be guided as the poten- 
tial depth is too small compared to their kinetic energy. 
Higher frequencies allow stable trajectories for molecules 
with higher Stark shifts and so the number of guidable 
states increases. Furthermore the guiding efficiency for 
molecules with higher Stark shifts is much higher. Both 
effects cause the rise in the count rate for frequencies 
below 8 kHz. For higher frequencies the count rate de- 



creases because the guided flux for every molecular state 
depends on the depth of the potential which drops off 
with uo~ 2 . 

Our conclusions are supported by two additional mea- 
surements. First, the dependence between the applied 
voltage and the frequency coo for which optimum guiding 
exists is analyzed. In the harmonic approximation, the 
stable region is given by |OT| < 3.75 with ft oc \ffi. As 
the curvature (3 scales linearly with the applied voltage 
V, we obtain ujq oc y/V. Experimental data are shown 
together with the expected square root dependence in the 
inset of Fig. 0Jl Unfortunately, the measurements cover 
only a small voltage interval, because for small voltages 
no guiding signal could be obtained, whereas sparks oc- 
curred for higher voltages. Nevertheless, the measured 
data agree with the calculated scaling. This allows one 
to choose the optimum switching frequency for a given 
voltage. In a second measurement, the flux is measured 
as a function of the applied voltage, using the optimum 
frequency for each voltage. The experimental data in 
Fig. ^3 show a quadratic rise in flux with increasing volt- 
age. This is expected for ND3, where most states have a 
linear Stark shift 

To summarize, two-dimensional trapping of neutral 
molecules with alternating electric fields is demonstrated 
experimentally. Our result proves that molecules with 
transverse and longitudinal temperatures of 20 mK and 
0.5 K, respectively, can be filtered out of a room tempera- 
ture effusive source efficiently. Our simulation shows that 
both HFS and LFS are guided, whereas the low-field seek- 
ers are guided and detected with a higher efficiency. The 
guiding efficiency for both classes will further increase if 
the injected molecules are precooled with cryogenic meth- 
ods. An important aspect of our guide is that it is suited 
to trap laser-cooled atoms as well 11]. Two-dimensional 
trapping of both atoms and molecules should therefore 
be possible. Finally, an extension of our technique to 
trap sufficiently slow molecules in three dimensions 
should be feasible. 
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